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Abstract. Consider a N x n matrix Z n = (^j\j 2 ) where the individual entries are a 
realization of a properly rescaled stationary gaussian random field: 

z ?ii2 = -7= E Hki,k 2 )U(h - fci.jfa - fca), 
(fci,fc 2 )e« 2 

where h a I 1 (<^ 2 ) is a deterministic complex summable sequence and (U (ji , 32); (ji , J2) £ Z 
is a sequence of independent complex gaussian random variables with mean zero and unit 
variance. 

The purpose of this article is to study the limiting empirical distribution of the eigen- 
values of Gram random matrices such as £T n Z* and [Z n + A n )(Z n + A n )* where A n 
is a deterministic matrix with appropriate assumptions in the case where n — > 00 and 

£-ce(o,oo). 

The proof relies on related results for matrices with independent but not identically 
distributed entries and substantially differs from related works in the literature (Boutct 
de Monvel et al. 0, Girko 0, etc.). 
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1. Introduction 

The model. Let Z n = (-Z^,- a , < ji < N, < j 2 < n) be a, N x n random matrix with 
entries 

Z hh = -7= E h (ki,k 2 )U(j 1 - h,j 2 - fc 2 ), 

where {U(j\, j 2 ), (iijj2) € Z 2 ) is a sequence of independent complex Gaussian random 
variables (r.v.) such that EU(j 1 ,j 2 ) = 0, EU(j 1 ,j 2 f = and E \U(ji, j 2 )\ 2 = 1, and 
(h(kx, k 2 ), (fci, k 2 ) £ Z 2 ) is a deterministic complex sequence satisfying 

\ h ( k i,k 2 )\ < 00 . 

(fci,fc 2 )ez 2 

The bidimensional process Zj 1 j 2 is a stationary gaussian field. Indeed, cov(Z™ ij9 , Z™,^, ) = 
n- 1 C(ji - j[,j 2 - f 2 ) where 

C(ji,.n)= KkiM)h*{ki-jx,k 2 -h) (1.1) 

(fei,fe 2 )ez 2 
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(we denote by a* the complex conjugate of a G C - we also denote by A* the hermitian 
adjoint of matrix A). 

The main results. The purpose of this article is to establish the convergence of the empirical 
distribution of the eigenvalues of various Gram matrices based on Z n . More precisely, we 
shall study the convergence of the spectral distribution of Z n Z* and (Z n + A n )(Z n + A n )* 
where A n is a deterministic matrix with a given structure. In particular, if Z n is square, 
we take A n to be Toeplitz. The contribution of this article is to provide a new method to 
study Gram matrices based on Gaussian fields. The main idea is to approximate the matrix 
Z n by a matrix Z n unitarily congruent to a matrix with independent but not identically 
distributed entries. This method will allow us to revisit the centered case Z n Z*, already 
studied by Boutet de Monvel et al. in P] and to establish the limiting spectral distribution 
of the non-centered case (Z n + A n )(Z n + A n )* for some deterministic matrix A n . 

Motivations. The motivations for such a work are twofold. First of all, we believe that this 
line of proof is new. Let us briefly describe the three main elements of it. 

The first one is a periodization scheme popular in signal processing and described as 
follows: 

Zn = fe 2 ) where Z? h = -= ^ h{k u k 2 )U ((jj - k x ) mod N, (j 2 - k 2 ) mod n) , 

^ n (/c 1 ,fc 2 )GZ 2 

where mod denotes modulo. 

The second element is an inequality due to Bai £Q involving the Levy distance C between 
distribution functions: 

C 4 (F AA ' ,F BB *) < -^-TrM — B)(A — B)*Tr(AA* + BB*), 

where F AA denotes the empirical distribution function of the eigenvalues of the matrix AA* 
and Tr(X) denotes the trace of matrix X. With the help of this inequality, we shall prove 
that Z n and Z n have the same limiting spectral distribution. 

The third element comes from the advantage of considering Z n . In fact, Z n is congruent 
(via Fourier unitary transforms) to a random matrix with independent but not identically 
distributed entries. Therefore, we can (and will) rely on results established in |Hj for Gram 
matrices with independent but not identically distributed entries. 

The second motivation comes from the field of wireless communications. In a communica- 
tion system employing antenna arrays at the transmitter and at the receiver sides, random 
matrices extracted from Gaussian fields are often good models for representing the radio 
communication channel. In this course, the stationary model as considered above is often 
a realistic channel model. The computations of popular receiver performance indexes such 
as Signal to Interference plus Noise Ratio or Shannon channel capacity heavily rely on the 
knowledge of the limiting spectral distribution of matrices of the type Z n Z^ (see 0,^01 an d 
also the tutorial JT] for further references). 

About the literature. Various Gram matrices based on Gaussian fields have already been 
studied in the literature. The study of the general case (Z n + A n )(Z n + A n )* has been 
undertaken by Girko in jjj. Since no assumptions are done on the structure of A n , there 
might not be any limiting spectral distribution. Girko finds asymptotic approximations of 
the Stieltjes transform of [Z n + A n )(Z n + A n )* . The method developed in is based on an 
exhaustive study of each entry of the resolvent {{Z n + A n )(Z n + A n )* — zl) -1 added to the 
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property that sufficiently remote entries arc asymptotically independent. 

Boutet de Monvel et al. [3j have also studied Gram matrices based on stationary Gaussian 
fields in the case where the matrix has the form V n + Z n Z*, V n being a deterministic Tocplitz 
matrix. Their line of proof is based on a direct study of the resolvent, taking advantage of 
the gaussianity of the entries. 

Disclaimer. In this paper, we study in detail the case where the entries of matrix Z n arc 
complex. In the real case, the general framework of the proof works as well if one considers 
the real counterpart of the Fourier unitary transforms, however the computations are more 
involved. We provide some details in Sectional 



2. Assumptions and useful results 

2.1. Notations, Assumptions, Stieltjes transforms and Stieltjes kernels. Let N = 

N(n) be a sequence of integers such that 

hm ^ = c. 

n-^oo TL 

We denote by i the complex number \/— 1, by 1a (x) the indicator function over set A and 
by S XQ (x) the Dirac measure at point xq. A sum will be equivalently written as X)fc=i or 
Y2k=i~n' We denote by CAf(0, 1) the distribution of the Gaussian complex random variable 
U satisfying EU = 0, EU 2 = 0, and E \U\ 2 = 1 (equivalently, U = A + LB where A and B 
are real independent Gaussian r.v.'s with mean and standard deviation each). 

Assumption A-l. The entries {Z™ 1]2 , < ji < N, < ji < n , n > 1) of the N X n 

matrix Z n are random variables defined as: 

Z ln =-)= E Kki,k 2 )U{ 3l - h,j 2 - fei), 
v n (fc 1 ,fc 2 )ez2 

where {h{k\,k2), (fci,fc2) G Z 2 ) is a deterministic complex sequence satisfying 

fcmax = E \h{kl,kl)\<QO 

and (U(ji,j2), {31,32) G Z 2 ) is a sequence of independent random variables with distribution 
CW(0,1). 

Remark 2.1. Assumption (A-QJ is a bit more restrictive than the related assumption [3], 
which only relies on the summability of the covariance function of the stationary process. 

For every matrix A, we denote by F AA , the empirical distribution function of the eigen- 
values of A A*. Since we will study at the same time the limiting spectrum of the matrices 
Z n Z* (resp. {Z n + A n ){Z n + A n )*) and Z* n Z n (resp. {Z„ + A n )* (Z„ + A n )), we can assume 
without loss of generality that c < 1. We also assume for simplicity that N < n. 

When dealing with vectors, the norm || • || will denote the Euclidean norm. In the case of 
matrices, the norm || • || will refer to the spectral norm. Denote by C + the set C + = {z G 
C, Im(z) > 0} and by C{X) the set of bounded continuous functions over a given topological 
space X endowed with the supremum norm || ■ W^. 
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Let /j, be a probability measure over R. Its Stieltjes transform / is defined by: 

We list below the main properties of the Stieltjes transforms that will be needed in the 
sequel. 

Proposition 2.1. The following properties hold true: 

(1) Let f be the Stieltjes transform of \i, then 

- the function f is analytic over C + , 

- the function f satisfies: |/(z)| < j^-jy, 

- if z G C+ then f(z) G C+, 

- if fx(— oo, 0) = i/ien z G C + implies z f(z) G C + . 

(2) Conversely, let f be a function analytic over C + such that f(z) G C + if z £ C + 
and |/(z)||/m(z)| bounded on C + . //liniy—f+oo —iyf{iy) = 1, then f is the Stieltjes 
transform of a probability measure \x and the following inversion formula holds: 

I f b 

H([a,b])= lim - / Im/(£ + irj) d£, 

where a and b are continuity points of \x. If moreover zf(z) G C + if z G C + then, 

fi(Mr) = 0. 

(3) Let P„ and P be probability measures over K and denote by f n and f their Stieltjes 
transforms. Then 

„ v ™ 



Nz G C+, f n {z) ► f(z)) 



Denote by M.c{X) the set of complex measures over the topological set X. In the sequel, 
we will call Stieltjes kernel every application 

tt : C+ Mc(X) 

either denoted Tr(z,dx) or Tr z (dx) and satisfying: 

(1) Vz G C+, Mg G C{X). 



gdn z 



< 



Im(z) 

(2) Mg G C(X), J gdir z is analytic over C + , 

(3) Vz G C+, V.g G C(X) and g > then Im (/ gdn z ) > 0, 

(4) Vz G C+ V.g G C(Af) and g > then Im (z / g d?r z ) > 0. 

2.2. A quick review of the results for matrices with independent entries. In order 
to establish the convergence of the empirical distribution of the eigenvalues, we will rely on 
the results based on matrices with independent but not identically distributed entries. Let 
us recall here those of interest (the assumptions and the statements are based on 

Consider a N x n random matrix Y n where the entries are given by 

_ Hji/N,j 2 /n) 

3132 /r Jlh 
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where X£j 2 and <f> arc defined below. 

Assumption A-2. The complex random variables {X™^ ; < ji < N, < j2 < n , n > 1) 

are independent and identically distributed (i.i.d.). They are centered with~E\XJ i: j 2 \ 2 = 1 and 
there exists e > such that E|X™ J2 | +e < oo. 

Assumption A-3. The function : [0, 1] x [0, 1] — ► C is such that |$| 2 is continuous and 
therefore there exist a non-negative constant $ max such that 

V(x,y) g [0,1] 2 , 0< |$(x,y)| 2 <$L x <oo. (2.1) 

Theorem 2.2 (independent entries, the centered case [H]). Assume that (A\3jj and (A- 
0) hold. Then the empirical distribution of the eigenvalues of the matrix Y n Y* converges 
a.s. to a non-random probability measure fj, whose Stieltjes transform f is given by f(z) = 
I[o l] n z(d x )> where ir z is the unique Stieljes kernel with support included in [0, 1] and satis- 
fying 

V.geC([0,l]), fgdn,= f- g( ^ 3(u , f) jdu. (2.2) 

If one adds a deterministic pseudo-diagonal matrix A n to the matrix Y n , the limiting 
equation is modified and in fact becomes a system of equations. 

Assumption A-4. Let A n = (A^-) be a complex deterministic N x n matrix whose non- 
diagonal entries are zero. We assume moreover that there exists a probability measure 
H( du,dX) over the set [0, 1] xl with compact support 7i such that 



1 N 



1=1 



Denote by Ti. c the support of the image of probability measure H under the application 
(u, A) — > (cu, A) and by 1Z the support of the measure lr c>1 i(ckt) ® 6o(dX) where ® denotes 
the product of measure. The set H = H. c U 7?. will be of importance in the sequel (see also 
Remarks 2.4 and 2.5 in [S] for more information). 

Theorem 2.3 (independent entries, the non-centered case [Sj). Assume that (A\^, (A- 
0) and (A^ hold. Then the empirical distributions of the eigenvalues of matrices (Y n + 
A n )(F Tt + A n )* and (Y n + A n )*(Y n + A„) converge a.s. to non-random probability measures 
\x and fi whose Stieltjes transforms f and f are given by 



f(z) = / TT z (dx) and f(z) = / TX z {dx) 

JTi J H 

where tt z and % z are the unique Stieljes kernels with supports included in Ti. and Ti. and 
satisfying 

( gdir, = [ , H(du,dX) (2.4) 

J J -z(l + lm 2 (u,t)n(z,dt,dC))+ 1+cImHt X cu)7l{M0 
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f gdn z = c [ dic^X) H(du,d\) 

f 1 gM) . 

( C 'J C -z(l + cJ\t>mt,u)TT(.Z,dt,dO) U 

(2.5) 

where \2.$ and \2. 5)) hold for every g € C(H) 



3. The limiting distribution in the centered stationary case 
We first introduce the following complex-valued function $ : [0, 1] x [0, 1] — > C denned by: 

<t>(h,t 2 ) = J2 Ki u i 2 )e 2 * i(lltl - l2t2) (3.1) 

((i,i 2 )ez 2 

We also introduce the p x p Fourier matrix F p = {F? i j 2 )o<ji,j 2 <p defined by: 

K , h = -Le^(lf). ,3.2, 
Note that matrix F p is a unitary matrix. 

Theorem 3.1 (stationary entries, the centered case j3J |7|). Let Z n be a N x n matrix 
satisfying (A^Tjjj. Then the empirical distribution of the eigenvalues of the matrix Z n Z^ 
converges in probability to the non-random probability measure pi defined in Theorem \2."2\ 

3.1. Proof of Theorem [TTJ Recall that 

Z lh = \ E Kki,k 2 )U{ 3l - h,h - to). 

We introduce the N x n matrix Z n whose entries are defined by 

Zjih = "7= h ( k ii k i)U{ji - fci mod N,j 2 - k 2 mod n). 

^ n (fei,fe2)ez 2 

For simplicity, we shall write U n (ji,j 2 ) instead of U(ji mod iV, j 2 mod n). Recall that C 
stands for the Levy distance between distribution functions. The main interest in dealing 
with matrix Z n lies in the following two lemmas. 

Lemma 3.2. Consider the N x n matrix Y n — FnZ u F*. Then the entries YJ™ ofY n can 
be written 

where $ is defined in IS. 1)) and the complex random variables {Xf^ , < l\ < N, < I2 < n} 
are independent with distribution Cj\f(0, 1). 
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Proof of Lemma \y.<H We first compute the individual entries of matrix Y n — F^Z n F*: 

32 = : n — 1 

-. 2irr( il ^ 1 i2 ' 2 ) 

= -p E £ m= " E Kkuk2)U n {h-ki,h-k2) 

v j 1 = o:iv-i VJV " (fc!.fe 2 )ez 2 

32 = : n — 1 
1 e 2i7r(itiL-i2^.) 

v ji = : AT - 1 v mi = : N - 1 

j2 = : n — 1 m,2 — : n — 1 

x E Mil -m i + ^l-^i J2 - m2 + k 2 n) 

(fei,fe 2 )6Z 2 

E ^(mi.ma)- 



^/n V N n J ^ y/Nn 

v x ' mi = : N - 1 v 

m 2 = : TO — 1 

Let be the random variable defined as 



\/A/Y) 

mi=0:JV-l VJV " 

m 2 = : n — 1 

for < l\ < N — 1 and < I2 < n — 1. Denoting by X„ and J7„ the N X n ma- 
trices with entries ^;™j 2 and U(li,l?) respectively, we then have JT„ = F^U n F*. Define 
vec(A) to be the vector obtained by stacking the columns of matrix A. Then the Nn x 1 
vectors X = vec(X„) and U = vec(f7 n ) are related by the equation X = (F* ® Fn)U 
(Lemma 4.3.1 in [Hj), where ® denotes the Kroncckcr product of matrices. The vector 
X is a complex Gaussian random vector that satisfies EX = (F* <g> Fn) EU = and 
EXX T = (F* <g> F N ) EUU T (F* <g> Fjv) = 0. After noticing that the matrix (F* <g> Fjy) is 
unitary, we furthermore have EXX* = (F* ® F/v) EUU* (F* ® F n )* = I nN where I p is the 
p x p identity matrix. In short, the entries of X n are independent and have the distribution 
CAA(0, 1). Lemma O is proved. □ 

Lemma 3.3. Let B n be a N x n deterministic matrix such that the sequence iTr_B„B* is 
bounded. Then 

£fp(z n +B n )(z„+B n )*^ F (z n +B n )(z n +B n )*\ p > 0j 

\ / n — >oo 

p 

where — > denotes convergence in probability. 
Proof of Lemma HOI Bai's inequality yields: 

ci{F i.Z n+ B n ) ( Z n+ B n Y iF (Z n+ B n) (Z n+ B nY) < 2 ^ _ f _ f 

x Tr [{Z n + B n )(Z„ + B n )* + (Z n + B n )(Z n + B n )*) (3.3) 
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Wc introduce the following notations: 



OL n = —Tr(Z n — Z n )(Z n -Z n )*, 
n 

(in = -Tr(Z n + B n ){Z n + B n )*, n = -Tr(Z n + B n )(Z n + B n f . 
n n 

With these notations, Inequality l|3.3[l becomes: 

£ 4 r F (Z n+ B n )(Z n+ B n r ;F (Z n +B n )(Z n+ B n r^ < 2a n (f3 n +P n ). 

In order to prove that £(f( z ™+ s "H z ™+ s ")* , F^n+B n )(z n +B n y-^ Q; it ig su ffi cient to prove 
that a n (P n + (3 n ) — > 0, which follows from a n — ► and (3 n and f3„ being tight. Indeed, 



a n > -^1 +P{/3„ > 2K}. 
2K J 



P{a„(/3„ + /3„) > e} < P{a„/3„ > e/2} + F{a„p„ > e/2} 
< P{a„ > ^} +P{/3„ > 2X} + p| 

Let us first prove that 

On A 0. (3.4) 
Since a n is non-negative, it is sufficient by Markov's inequality to prove that Ea„ — > 0. 



*n — — Tr(Z„ — Z n )(Z n — Z nj 



= - y 

31 = : N - 1 

32 — : n — 1 

= - y 



771 _ hn 
31,32 ~ 31,32 



ji = : JV - 1 (fci,fc 2 
32 = : n — 1 



y h(k 1 ,k 2 )V(j 1 - ki,]2 - k 2 ) 



where V(ji,j 2 ) stands for U(ji,j 2 ) — U n (ji, j 2 ). Thus 

31 = 0: iV- 1 (fei,fc 2 ) 6 Z 2 

32 = : n- 1 j/) s z 2 

Introduce the set J" = {0, • • ■ , AT - 1} x {0, ■ • • , n - I}. Then 



iV(l 1 ,l 2 )V*(l' 1 ,l' 2 )= l Z 2_j(h,l 2 ) 1 Z 2_ j(l[,l' 2 ) 



+ + miN,l' 2 + m 2 n) 

(mi,m2)G^ 2 
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and Ea„ becomes Ea„ = Ea„j +Ea„, 2 where 

E\l = ~2 E X] l^( fc l> fc 2)| 2 Iz^-jO'l ~ ^1, J2 - ^2), 



n — 
n=0:iv-i (fc 1 ,fe 2 )ez 2 

32 = : n - 1 



Ea„,2 = — ^ X k 2)h*(k[, fc 2 ) !z 2 -j0'i ~ fci,j2 - h) 



n 

31 = : IV - 1 (fci, fc 2 ) e z 2 
J 2 = : n - 1 (k[, k' 2 ) 6 Z 2 



x lz 2 -j(h - Kih - k 2 ) X 1 (fei,fc2)( fe 'i +miAT,4 + m 2 n) 

(mi,m 2 )GZ 2 

Let us first deal with Ea„,2. 

Ea„,2 < ^ E X \K k uki)\ lz 2 -j(ji - fci, j 2 - k 2 ) 

J! = : JV - 1 (fcilsJeZ 2 
32 = : n — 1 

x J2 \ h (kiX)\^-j{h-k'i,h-k 2 ) 

x E 1 (fci,fc2)( fc i +rniiV, k' 2 + m 2 n). 

{ra\^ra 2 )^ r L' 1 

Since /i is summable over 1? by (A-QJ, 

E l^( fc i! fc 2)l ^-^(Ji-MjJ^-^) X 1 (kiM) (K +m 1 N,k 2 +m 2 n) 
(fei,fc^)ez 2 (m 1 ,m 2 )eP 

is bounded by /i max and 

EV<^ E E l ft (*i.*a)l lp-jUi-ki,j 2 -k 2 ). (3.5) 

Since 



n z 

31 = : at - 1 (&i,fc 2 )eJ 

32 = : n - 1 



we get: 



1 / • , ■ 7 \ / ji - fci < or j'i - ^ > JV, 

^ |/i(fci,fc 2 )| lz 2 -j(ji - fci, h - k 2 ) 





(fei,fe 2 )ez 2 








E i^(fci,fc 2 )i + 




E IMfci,fc 2 )l 




fc^ — —00 : j\ — N; 
k 2 — — 00 : j 2 ~ n 


fci = 
fe 2 = 


— 00 : — JV; 
J2 + 1 : ao 




+ IMfci,fc 2 )H 




E IMfci,fc 2 )|. 




fci = 31 + 1 : 00; 
k 2 — —00 : j 2 — n 


fci 
fc2 


= 31 + 1 : 00; 
= j '2 + 1 : oa 


The changes 


of variable | ^} ^ Jl and < 


r J2 
1 fc 2 


= n — 1 — 7o . , , 

= -fc 2 yldd 


E 


E IMfci,fc 2 )l = 


E 


E IM-fci,-fc 2 


31 = : N 

32 = : n - 


- 1 fci = 00 : 31 - JV; 3' = 

- 1 fc 2 = — 00 : 32 — n 32 - 


: : AT 
: : n 


— 1 fcj_=ji + l:ooj 

- 1 fc 2 — 32 + 1 : 00 
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By performing similar changes of variables, one gets: 



E E \ h ( k i,k 2 )\ lz 2 -j(ji ~fa,h ~ k 2 ) 

J1=0:JV-1 (fci,fc 2 )GZ 2 
32 = : n — 1 

E E \H-k 1 ,-k2)\ + \h(-k 1 ,k 2 )\ + \h(k 1 ,-k 2 )\ + \h(k u k 2 )\ 



31 — : JV — 1 fci — 3*i + 1 : oo; 
j 2 = : n — 1 /C2 — J2 + 1 : oo 



S(J1J2) 

In order to check that 



- 2 E ^'i'^) ►O, (3.6) 

3! = : JV - 1 
32 = : n — 1 

we introduce T(j) = Efc 1 +fe 2 > J +2 l^( _fc x, ~k 2 )\ + \h(-ki, k 2 )\ + \h(ki,-k 2 )\ + \h(k 1: k 2 )\. Is 
is straightforward to check that T(j) > and that S(ji,j 2 ) < T(j\ + j 2 ). We prove 

(|3.6|) by a Cesaro-like argument: Let no < N be such that T(uq + 1) < e. Wc have 

\ E S{ji 1 j 2 ) = \ V S(ji,j 2 ) + ~2 V S(j lt j 2 ). 

31 = : N - 1 0<ji+j 2 <no n + 1 < 31 + 32\ 

32 — : n — 1 31^-^™l)32^; n — 1 

(3.7) 

If n is large enough, then the first part of the right handsidc of 13.7|l is lower than e. Moreover, 

^ E s (n,h)<^ E n«o+i)<e 

"0 + 1 < Ji + J2 ; »o + i < 31 + 32 ; 

31 < N — 1, 32 < n— 1 ji < N — 1, 32 < n — 1 

and H3.6fl is proved. By pluging 1)3. 6[) into (|3.5|l . we prove that Ea„,2 — * 0. Using the same 

kind of arguments, one proves that Ea„.i — > 0. Finally, l|3.4|l is proved: a n — > 0. 
Let us now check that 

3K>0, E/3 n <K and 3K > 0, Ep n <K. (3.8) 

This will imply the tightness of (i n and (3 n . 

Recall that by assumption there exists -B max such that sup„ — Tr B n B^ < B max . Consider 
now: 

(i i\ 2 

(-TiZ n Z*Y + (-TiB n B* 
\n J \n 

In particular, 

Tr ( Z n + B n )( Z n + B n Y TrZ n Z* I „ ( TrZ n Z* \ * 

n n \ n J 

(°) „TvZ n Z* I / {TrZ n Z*\\^ , x 

< E ^ + 2B£ax(E[ + #max (3.9) 



LIMITING DISTRIBUTION OF STATIONARY GRAM MATRICES 



11 



where (a) follows from Jensen's inequality. Notice that (|3.9|) still holds if one replaces Z n by 
Z„. Therefore in order to prove (|3.8|) . it is sufficient to prove that: 



3 K' > 0, E 



( TrZ n Z* 



< K' and 3 K' > 0, E 



Tr Z n Z* 



< K . 



Consider 



E 



/Ti-7 7*\ 1 AT 

^ = - y, Ei^ J2 i 2 = vviE^rii 2 = -c(o,o), 



V n 



j 1 = l:N 
32 = 1 : n 



where C is defined by p. 1(1 . This quantity is asymptotically bounded. From lemma l3~2l we 
have 



E 



Tr Z n Z^ 



E 



Try y* 



1 



E 



31 = 1 : JV 
j 2 = 1 : n 



N' n 



\X n . I 2 < — $ 2 

' 3132 I — n ^mi 



which is also asymptotically bounded. Eq. (|3.8|) is proved and so is Lemma 13.31 



□ 



Proof of Theorem Vd.lX Lemma 13.31 implies that 

P {£ (> z " z ™, F^'A > e| ► for every e > 0. (3.10) 

By lemma T3.2I FpfZ n Z*F£f = Y n Y*. Since Fn is unitary, Z n Z* and K^^t have the same 
eigenvalues. Moreover, matrix Y n fulfills (A-J2J) and the variance profile $ defined in (|3.1|l 
satisfies (A-|5J since (/i(fci,&2) G (fcijfe)^ 2 ) is summable; therefore one can apply Theorem 
12.21 In particular, 

F Zr,z' n > ^ a. s . Ve>0, P {/I (f z " z - , ^ > e} >0 (3.11) 

where [i is the probability distribution defined in Theorem 12.21 Eq. 13.10f) together with 
(|3.11|) imply that F Zr,z ^ fi and Theorem 13. II is proved. □ 



4. The limiting distribution in the non-centered stationary case 

Recall the definitions of function <!> and matrix F p (respectively defined in l|3.1|) and 13.2|l ). 

Theorem 4.1 (stationary entries, the non-centered case). Let Z n be a N x n matrix 
satisfying fA™): let A n be a N x n matrix such that A„ = F]yA n F* is N X n pseudo- 
diagonal and satisfies (A^. Then the empirical distributions of the eigenvalues of matrices 
(Z n + A n )(Z n + A n )* and (Z n + A n )*(Z n + A n ) converge in probability to the non-random 
probability measures [i and fi defined in Theorem \2.cA 

Proof of Theorem \4~l\ We denote by F n = F^+ A "^ z "+ A "'>" and F n = F^ z "+ A '^ z " +A "'>" . 
Since A„ satisfies (A-0J, iTrA n A* = iTrA n A* is bounded and Lemma f3 . 31 implies that 

P(|£(F",F")| > e\ >0 for every e > 0. (4.1) 

By lemma HOI and the assumption over A n , 

(Z n + A n )(Z n + A n )* = F N (Y n + A„)(F„ + A n )*F* N . 
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Since the Fourier matrix is unitary, [Z n + A n )(Z n + A n )* and (Y n + A„)(l^, + A ra )* have 
the same eigenvalues. Since $ defined in i|3.1|) satisfies (A-|3J, the matrices Y n and A„ fulfill 
assumptions (A-|2J), (A-|3J) and (A-QJ therefore one can apply Theorem 12.31 In particular, 

F n ► n a.s. => Ve>0, P (\C(F n , fi)\ > e) >0 (4.2) 

where fi is the probability distribution defined in Theorem 12. 31 Eq. <|4.1|) together with (|4.2|l 
imply that F n — > P and Theorem 14. II is proved. □ 

In the square case n x n, we can deal with slightly more general matrices A n . 

Assumption A- 5. The n x n matrix A n is a Toeplitz matrix defined as A n = (a(ji — 
j2))o<ji,j2<fi w here (a(j'))jgz is a deterministic sequence of complex numbers satisfying: 

J2wj)\ <°°- 

iez 

Let tp : [0, 1] i—* C be the so called symbol of A n defined as 

</>(*) = $>(i)e 2i ^. (4.3) 

Due to (A-EJ, ip is bounded and continuous. 

Theorem 4.2 (stationary entries, the non-centered square case). Let Z n be a n x n matrix 
satisfying (A^Djl: let A n be a n x n matrix satisfying fAWjl. Then the empirical distributions 
of the eigenvalues of matrices (Z n + A n )(Z n + A n )* and (Z n + A n )*(Z n + A n ) converge in 
probability to non-random probability measures fi and ft whose Stieltjes transforms f and f 
are given by 

f(z) = / ir z (dx) and f{z)= / Tr z (dx) 

J [0,1] J [0,1] 

where n z and tt z are the unique Stieltjes kernels with supports included in [0, 1] and satisfying 
the system of equations: 

9d** = I ~ . — . W du ( 4 - 4 ) 



o -^(i+jw«,oi a ^)+ 1+J CTi^. 

q 6fk z = / 1 i , . ,, du (4.5) 



for every function g € C([0, 1]). 

Proof. The proof is based on the fact that a Toeplitz matrix A n is very close to a Toeplitz 
circulant matrix A n defined in such a way that the diagonal matrix A n = F n A n F* satisfies 
assumption (A-4). Denoting by ipn the truncated function ij) n (t) = Y^=- n a U) e 
choose A n to be the matrix whose entries are defined by 

n-1 



k=0 



Notice that in this case, A„ = F n A n F* is given by A„ = diag ( [i/j n (0), ip n (±),..., V'n( n ^ i )] ) 
where diag(w) is the diagonal matrix bearing the entries of the vector v on its diagonal. 
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One can also prove that the complex number a n {j\ — j'2) = ® 1 j l j 2 satisfies a n (0) = a(0) 
a(n) + a(— n) and 



S B (?) 



+ a ti -n) if n - 1 > j > 0, 
+ a 0' + if - n + 1 < j < 0. 



We denote by F n and F" the distribution functions F n = F^ z ^ +A ^ z ^ +A ^* and F n = 
F (z n +A n )(z n +A n )\ We ghall proye that c{F n 1 F n ) -► as n -> 00. 

Bai's inequality yields: 

£ 4 (F", F n ) < -^(A, - in)(A, - i„)*Tr(A n ^ + i„I*). (4.6) 

We first prove that n~ 1 Ti:(A n A* n ) and n~ lr £i(A n A* n ) are bounded: 

i 1 n ~ * n— * / I M\ 2 

- t ^«a; = - E KA-i 2 )i 2 = E Ia(3)l a (l-^)<(EI°CJ)|) ■ ( 4 -7) 

iij2=o j=— n+i v / jez 

Moreover, 



n-l 



n n n 



-TtA n A* n = — TrA„A* = - E M - ) < ( E ( 4 ' 8 ) 



2 



We now prove that 

-Tr{A n -A n ){A n -A n )* ► 0. (4.9) 

n — >oo 

Indeed, 



1 1 

-lr(A„-I n )(A n -!„)* = - E l«0'i-J2)-5"(j 1 -j 2 )| 5 

3l ,32=0 

E Wi)-a"0-)l 2 (i-^) 



j=-(n-l) 

n-l 



n—x / 

a(-n) + a(n)\ 2 + E (Hj ~ n)| 2 + |o(n - j)| 2 ) [l 
a (-n)+a(n)| 2 + E-(|a(j)| 2 + l«(-j)| 2 



"i)| 2 



n-l 

n 

-, J oo 

< | a (_n)+ a( n)| 2 + -^ J (|a( J )| 2 + |a(^)| 2 )+ E (V(i)| 2 + K 

3 = 1 j=J+l 

By first taking J large enough then n large enough, the claim is proved by a 2e-argument. 
Eq. I)4.6|l together with the arguments provided by (|4.7(l . 1)4. 8[) and 1)4. 9J) imply that 

C(F n ,F n ) >0. 

n — >oo 

It remains to prove that F" converges towards the non random probability distribution 
characterized by equations (|4.4(l and (|4.5|l . As previously, the variance profile $ defined in 
(13. 1|) satisfies (A-|3J). Moreover, we have 



1 ™ 

i—1 
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where H(du,dX) is the image of the Lebesgue measure over [0,1] under u h- > (u, \ip(u)\ 2 ). 
Therefore A n satisfies (A-4) and Theorem 14. II can be applied. This completes the proof of 
Theorem O □ 



5. Remarks on the real case 
In the case where the entries of matrix Z n are given by 

^ n (fci,fc2)ez 2 

where (/i(fci,fc2), (^1,^2) G Z 2 ) is a deterministic real and summable sequence and where 
U (31,32) are real standard independent gaussian r.v.'s, the conclusion of Lemma |3.2I is no 
longer valid. In fact the entries of Y n = FnZ u F* are far from being independent since 
straightforward computation yields: 



I N-U ,n-l 



for < h < N and < l 2 < n. 
We introduce the p x p orthogonal matrix Q p = {Qj' 1 j 2 )o<j 1 ,j 2 <p defined by: 



) p - 



1 



VP' 



0<h< p- 



In the case where p is even, the entries Q p (ji, 32) (ii > 1) are defined by 

fl cos (^f^ if l<ji<f-l > 0<j 3 <p; 
/|sin( 2 ^) if 1< ji < f-1,0 <j 2 <p; 



» 231-1,32 

) P = 

5 231,32 

= (-!) J 

J2 



Q^-i 



if < j a < p. 



In the case where p is odd, they are defined by 




sin 



os(^™) if l<ji<^,0<j 2 <p; 
j^) if l<ii<^,0<i 2 <p. 



In the sequel, [^J stands for the integer part of x. The following result is the counterpart of 
Lemma 13.21 in the real case. 

Lemma 5.1. Consider the N x n matrix W n = Qjv^nQn where A T is the transpose of 
matrix A. Then the entries Wfo of W n can be written as 



1 


*(- 


h + l 


1 


h + l 




\fn 


2 


n 


2 


) 



V7J 



where $ is defined in \S.1\) and the real random variables ,0 < l\ < N, < l 2 < n} 

are independent standard gaussian r.v. 's. 

The proof is computationally more involved but similar in spirit to that of Lemma 13.21 It 
is thus ommited. 

As a consequence of this lemma, Theorems 13.11 and 14.11 remain true with the following 
minor modification: In Eq. H2.2JI. I|2.4|l and (|2.5|l . the quantity |3>| 2 must be replaced by $^ 
where 

- |*(«/2,«/2)|. 
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Similarly, in the case where the Toeplitz matrix A n introduced in (A-[SJ) is real, Theorem 14. 21 
remains true if one replaces in 1)4. 4|l and (|4.5|) the quantities |$| 2 and \ip\ 2 by and t/'r 
where 

M«) = l^(«/2)|. 

The proof of Theorem 14.21 can be modified by replacing the Fourier matrices F p by Q p (see 
also 0], chap. 4 for elements about the pseudo-diagonalization of a real Toeplitz matrix via 
real orthogonal matrices Q p ). 
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